Weakly bound dimers of fermionic atoms 
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We discuss the behavior of weakly bound bosonic dimers formed in a cold Fermi gas at a large 
positive scattering length a for the interspecies interaction. We find the exact solution for the dimer- 
dimer elastic scattering and obtain a strong decrease of their collisional relaxation and decay with 
increasing a. The large ratio of the elastic to inelastic rate is promising for achieving Bose-Einstein 
condensation of the dimers and cooling the condensed gas to very low temperatures. 
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The last years were marked by remarkable achieve- 
ments in the physics of cold Fermi gases. Several groups 
succeeded in cooling trapped fermionic atoms to well be- 
low the temperature of quantum degeneracy (Fermi tem- 
perature Tp) [1-7], and the ratio T/Tp — 0.05 has been 
reached in the recent MIT experiment [8]. One of the 
main goals of these studies is achieving a transition to 
a superfluid Cooper-paired state. Trapped Fermi gases 
are very cold and dilute, with temperatures T < 1 /iK 
and densities n ~ 10 13 cm -3 , and for an attractive in- 
terspecies interaction (negative s-wave scattering length 
a) the most efficient will be the superfluid s-wave pairing 
between atoms of different components. However, the su- 
perfluid transition temperature T c is exponentially small 
compared to the Fermi temperature Tp and is beyond 
experimental reach for ordinary small values of a. 

At present, actively discussed ideas to circumvent this 
difficulty rely on superfluid pairing between fermionic 
atoms via a Feshbach resonance [9-11]. In the vicin- 
ity of the resonance the scattering length is very large, 
being negative on one side of the resonance and positive 
on the other side. On approach to the resonance, the 
gas enters a strong-coupling regime. This occurs when 
the parameter kp\a\ > 1, with kp = \ / 2mTp/h being 
the Fermi momentum and m the atom mass. Crossing 
the resonance and making the scattering length positive, 
the formation of weakly bound dimers of two different 
fermions becomes energetically favorable. Sufficiently far 
from resonance on the positive side, one has a weakly in- 
teracting gas of these composite bosons and encounters 
the problem of their Bose-Einstein condensation (BEC). 
This cross-over to the BEC regime has been discussed in 
literature in the context of superconductivity [12-14]. 

The BEC regime of the bosonic dimers is interesting 
from a fundamental point of view as it couples the prob- 
lem of superfluidity in Fermi gases to the problem of 
molecular condensates. Here the most important ques- 
tions are the stability of the condensate with regard to 
elastic dimer-dimer interactions, and the decay of the gas 
due to relaxation of the dimers to deep bound states. The 



relaxation occurs in dimer-dimer collisions and in colli- 
sions of dimers with remaining unbound atoms. It is a 
crucial process as these dimers are diatomic molecules in 
the highest rovibrational state. Several experiments show 
that such molecules consisting of bosonic 87 Rb [15,16] 
and 133 Cs atoms [17], or fermionic 40 K atoms with a 
scattering length a ~ 100A [18], undergo a rapid col- 
lisional decay. On the other hand, recent observations 
[19-21] indicate the existence of long-lived weakly bound 
LI2 dimers at densities ~ 10 13 cm -3 . 

In this Letter we present an exact solution for the 
dimer-dimer elastic scattering, assuming that the (pos- 
itive) scattering length greatly exceeds the characteristic 
radius of interaction between atoms, R e . Then, as in the 
case of the 3-body problem with fermions (see [22,23] and 
references in [23]), the amplitude of elastic interaction is 
determined by a single parameter, the two-body scatter- 
ing length a, and can be found in the zero-range approx- 
imation for the interatomic potential. Our findings lead 
to a positive dimer-dimer scattering length add — 0.6a 
[24] . This is quite different from the assumption of ear- 
lier studies, add = 2a [14], and has serious consequences 
for the low-temperature behavior of the system. 

We then discuss the collisional relaxation of the weakly 
bound dimers to deep bound states and show that it is 
suppressed due to Fermi statistics for the atoms. The 
binding energy of the dimers is £0 = h 2 /ma 2 and their 
size is close to a. The size of deep bound states is of 
the order of R e <C a. Hence, the relaxation requires 
the presence of at least three fermions at distances ~ R e 
from each other. As two of them are necessarily identical, 
due to the Pauli exclusion principle the relaxation prob- 
ability acquires a small factor proportional to a power 
of (fci? e ), where k ~ 1/a is a characteristic momentum 
of the atoms in the weakly bound molecular state. The 
inequality a^> R e allows us to obtain the dependence of 
the atom-dimer and dimer-dimer wavefunctions at short 
interparticle distances on the two-body scattering length 
a and thus to establish a strong decrease of the relax- 
ation rate with increasing a. Our results show that the 
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rate constant is proportional to a~ s , where s is close to 
3. This indicates a remarkable collisional stability of the 
weakly bound bosonic dimers at a large positive a, which 
is consistent with recent observations at ENS [19], Inns- 
bruck [20], and Rice [21]. Note that the elastic collisional 
rate is proportional to a 2 , and for large a we obtain a very 
large ratio of this rate to the inelastic rate of relaxation. 

The key idea behind these studies is twofold. First, a 
large ratio of the elastic to inelastic rate is promising for 
evaporative cooling of the dimers, achieving their BEC, 
and further cooling of the Bose-condensed gas to very low 
temperatures. Molecular BEC is a great challenge and 
many ongoing studies are directed to reaching this goal. 
Second, adiabatically crossing the Feshbach resonance 
and reaching a negative value of the scattering length 
a, the molecules are converted into fcrmionic atoms at 
extremely low temperatures T <~ 0.01TV which can be 
below the BCS transition temperature T c [25]. Remark- 
ably, at these temperatures the strong Pauli blocking of 
elastic collisions provides the collisionlcss regime for the 
atomic Fermi gas even for a large negative a. This opens 
up the possibility of identifying the BCS-paired state via 
the asymmetric expansion of the cloud released from a 
cylindrical trap [26] . 

We start with the dimer-dimer elastic scattering which 
is a 4-body problem described by the Schrodinger equa- 
tion 
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mE/h 2 ] # = -(m/h 2 ) 
U(n) + U(r 2 ) + J2 U (( r i + r 2 ± V2R)/2 



*■ (1) 



Here ri is the distance between two given distinguishable 
fermions, r 2 is the distance between the other two, R/ \pl 
is the distance between the centers of mass of these pairs, 
and U is the interatomic potential. The total energy is 
E = — 2e +£, with e being the collision energy. 

The wavefunction \P is symmetric with respect to the 
permutation of the dimers and antisymmetric with re- 
spect to permutations of identical fermions: 

*(n,r 2 ,R) = *(r 2 ,n,-R) (2) 
'r!+r 2 ±V2R n + r 2 T \/2R ±(n-r 2 )\ 
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At energies e <C £o the scattering is dominated by the 
contribution of the s-wave channel and can be analyzed 
from the solution of Eq. (1) with E = —2s . For large R 
the corresponding wavefunction is given by 

*KKpo{ri)Mr2)(l-y/2aM/R), (3) 
where add is the dimer-dimer scattering length, and 

<Po(r) = (r\ / 2ira)~ 1 exp (— r/a) (4) 
is the wavefunction of a weakly bound dimer. 



The characteristic de Broglie wavelength of atoms is 
of the order of a and it greatly exceeds the radius R e 
of the interatomic potential U. Hence, for finding the 
scattering amplitude one can use the well known method 
of pseudopotentials, that is, to replace the potential U by 
a pseudopotential providing proper boundary conditions 
for ^ at vanishing distances between two distinguishable 
fermions. For ri — > the boundary condition reads: 



*(ri,r 2 ,R) - /(r 2 ,R)(l/47rn - 1/4™), 



(5) 



where the function /(r 2 , R) carries the information about 
the second pair of particles when the first two sit on top of 
each other. The boundary conditions for r 2 — ► 0, and for 
ri+r 2 ±\/2R — * arc readily written using Eqs. (2). The 
knowledge of the function / is sufficient for describing 
all of the four boundaries. Combining Eqs. (3) and (5) 
one can deduce the value of add from the behavior of / 
at large R. It is important that in contrast to the 
function / depends only on three variables: the absolute 
values of r 2 and R, and the angle between them. In the 
following, we derive and solve the equation for /. 

In the pseudopotential approach the interaction po- 
tential can be written as U(r) — (4iTh 2 a/m)5(r)(d/dr)r. 
Then, making use of Eq.(5) we express the rhs of 
Eq.(l) in terms of the function /. The first term reads 
— mU(ri)'&/K — <5(ri)/(r 2 , R), and similarly for the 
other terms. The resulting inhomogeneous Poisson equa- 
tion can be solved by inverting the differential oper- 
ator. The Green function of this operator is unique 
(E = — 2e < 0), and is expressed through / as: 

*(S)= [[ ^ G(\S - Si\) G(\S-S ± \)] f(r', R'), (6) 

r / R ,i=l,2 ± 

where S = {n,r 2 ,R}, Si = {0, r', R'}, S 2 = {r\ 0,-R'}, 
S ± ={r'/2±R'/\/2,r'/2=F R'/V^Tr'/v^}, and the Green 
function is G(x) = {2ir)- 9 / 2 (xa/V2)- 7 / 2 K 7/2 (V2x/a), 
with Kf/2 being the decaying Bcssel function. 

The equation for the function / is obtained from 
Eq. (6) by taking the limit n — * and using Eq. (5). 
The singular terms proportional to 1 /n cancel each other 
automatically, and matching the regular parts yields 

J {G(\Si 5!|)[/(r', R') - /(r, R)] + [G(\Si - S 2 |) 

3/ 

- ^Gd^-SiD] /(r', R')} = (V2- l)/(r, R)/4™, (7) 

where Si — {0,r,R}. Note that the only length scale 
in Eq. (7) is a. Solving this equation numerically we 
obtain the function / for all distances R and r. Ac- 
cording to Eqs. (3), (4) and (5), at R — > oo one has 
/ = (2/ra) exp(— r/a)(l — \/2add/R)- Fitting our results 
with this expression we find with 2% accuracy: 
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a dd = 0.6a > 0. (8) 

Our calculations show the absence of 4-body weakly 
bound states, and the behavior of / at small R suggests a 
soft-core repulsion between molecules, with a range ~ a. 
The result of Eq. (8) indicates the stability of molecu- 
lar BEC with respect to collapse. Compared to earlier 
studies which were assuming add = 2a [14], Eq. (8) gives 
almost twice as small a sound velocity of the molecular 
condensate and a rate of elastic collisions smaller by an 
order of magnitude. 

The lifetime of the Bose gas of weakly bound dimers is 
determined by the process of their collisional relaxation 
into deep bound states. The released binding energy of 
a deep state is of the order of ft 2 /mR 2 . It is transformed 
into the kinetic energy of particles in the outgoing colli- 
sional channel and they escape from the trapped sample. 
We will establish a dependence of the relaxation loss rate 
on the scattering length a, without going into a detailed 
analysis of the short-range behavior of the systems of 
three and four atoms. We employ a perturbative scheme 
assuming that the amplitude of the inelastic process of 
relaxation is much smaller than the amplitude of elas- 
tic scattering. Then the dependence of the relaxation 
rate on the scattering length a is related only to the a- 
dependence of the initial-state wavefunction 'P. 

We first discuss the relaxation of weakly bound dimers 
to a deep bound state in their collisions with atoms. This 
process occurs when all of the three atoms approach each 
other to distances ~ R e . As at all interparticle distances 
r <C a the three-body wavefunction 'J is determined by 
the Schrodinger equation with zero collision energy, it 
depends on a only through a normalization coefficient: 
\I> = A(a)ip. In the region where R e <C r <C a, the o- 
independent function ip can be found in the zero-range 
approximation. Then the coefficient A(a), which deter- 
mines the dependence of the relaxation rate on a, is ob- 
tained by matching <J> with the solution at large interpar- 
ticle distances. 

This requires us to find <J> at distances ~ a, which we 
do using the zero-range approximation described in Ref. 
[23]. As in the case of four particles, we introduce the 
corresponding function /(r). The correctly symmetrized 
wavefunction \I>(x, y) is expressed through / as: 

* = E/ ±(5 (\A x - r / 2 ) 2 + A/3r/2)^ /», (9) 

r 

where y is the distance between identical fermions, 
\/3x/2 is the distance between their center of mass 
and the third atom, and the Green function G{z) = 
(8ir 2 z 2 a 2 )~ 1 K2(z/a). At a vanishing distance between 
two given distinguishable fermions r± = (\/3x± y)/2 — > 
0, the function \P (9) should have a correct asymptotic 
behavior in analogy with Eq. (5). This gives an equation 
for the function /(r): 



J G(k-r'l) [/>' ) - /> )] - G (sjv 2 + r ' 2 + rr'j />')= 0. ( 10) 

r' 

For small separations between all three fermions, 
Eqs. (10) and (9) give the same result as the method 
of hyperspherical harmonics [28]: W w A(a)p 7 <I> 7 (£l), 
where 17 is a set of hyperangles and the hyperradius 
p = \J 'x 2 + y 2 satisfies the inequality p <C a. The pa- 
rameter 7 depends on the symmetry of the wavefunc- 
tion and on the orbital angular momentum of the atom- 
dimer motion. For the s-wave atom-dimer scattering we 
have /(r) = /(|r|) and obtain 7 = 71 w 0.1662. At large 
atom-dimer distances y we have ^> w ±</?o( r ±)(l~ a ad/y), 
where ipg is given by Eq. (4), and a a d = 1.2a is the atom- 
dimer scattering length [27,23]. Solving Eqs. (10) and (9) 
at distances of the order of a and thus matching the two 
asymptotic solutions, we find A(a) = const xa~ 3 / 2 ~ 71 . 
This result is explained extending the wavefunction * to 
interparticle distances ~ a from the region of small dis- 
tances. This gives ^^A(a)a 71 , whereas being extended 
to distances ~ a from large separations, the wavefunction 
is ~ 1/ a 3 / 2 . Equalizing the two expressions leads to the 
obtained coefficient A(a). The numerical constant is not 
needed as the absolute value of the relaxation rate is de- 
termined by the contribution of interparticle distances 
~ R e , where the zero-range approximation is not valid. 
It only gives a correct dependence of the rate on a. 

As the relaxation rate constant a rc i oc A 2 (a) , we obtain 
a rc i oc a~ 3 ~ 271 = a~ 3 ' 33 . One thus sees that the relax- 
ation rate due to dimer-atom collisions rapidly decreases 
with increasing two-body scattering length a. 

The obtained results can be easily generalized to the 
case of the s-wave dimer-dimer scattering. Indeed, the 
relaxation process requires only three atoms to approach 
each other to short distances. The fourth particle does 
not participate. Let p and O be the hyperradius and 
the set of hyperangles of the system of three fermions 
and x the distance between their center of mass and the 
fourth atom. At p <~ R e and |x| ^> R e the four-particle 
wavefunction decomposes into p 72 $ 72 (J7)r/(x), where the 
function r/(x) describes the motion of the fourth particle. 
Averaging the relaxation probability over the motion of 
the fourth particle makes the problem similar to the re- 
laxation in atom-dimer collisions, and we arrive at the 
same result for the rate constant. However, there is a 
relaxation channel that is more important in the limit 
of very large a. For the s-wave dimer-dimer scattering, 
both the fourth particle and the atom bound to this par- 
ticle can undergo p-w&ve scattering on the other dimer in 
such a way that their total angular momentum is equal 
to zero. This corresponds to an antisymmetric function 
/(r) in Eqs. (9) and (10) and leads to 72 = -0.2273. For 
the relaxation rate we then find a ro i oc a~ 3 ~ 272 = a -2 ' 55 . 

Assuming that the short-range physics is characterized 
by the length scale R e and the energy scale ft 2 /mR 2 , 
we can restore the dimensions and write: a rc i = 
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C(TiR e /m)(R e /a) 8 , where s w 2.55 for the dimer-dimer 
collisional relaxation, and s 3.33 for the relaxation in 
atom-dimer collisions. The coefficient C depends on a 
particular system and is different for dimer-dimer and 
atom-dimer collisions. Obviously, in the limit R e /a — > 
the dimer-dimer relaxation should dominate over the 
atom-dimer one. 

The much slower collisional decay of the weakly bound 
dimers at larger a, following from our results, is consis- 
tent with recent observations for Li2 molecules [19-21], 
although the experimental data do not yet provide the 
dependence of a re i on a. Nevertheless, the ENS result 
a IC \ «2x 10 -13 cm 3 /s for the L12 loss rate at a « 1500ao 
[19], allows us to establish the coefficient C w 20. 

We emphasize that the remarkable stability of such 
weakly bound dimers at a 3> R e is due to Fermi statistics. 
Indeed, two identical fermions participating in the relax- 
ation process have very small relative momenta k ~ l/o 
and, hence, the process is suppressed compared to the 
case of a dimerized gas of bosonic atoms. 

The inelastic rate constant a ro i is much smaller than 
the rate constant of elastic collisions %na 2 dd VT, where vt 
is the thermal velocity. For the Li2 weakly bound dimers 
at a temperature T <~ 3fiK and a ~ 1500ao, the corre- 
sponding ratio is of the order of 10~ 4 or 10~ 5 . This is 
nearly the same as in a cold gas of Rb atoms at densities 
n <~ 10 14 cm~ 3 . For Li2 the values of both rate constants 
are larger, and one can expect even a higher efficiency 
of evaporative cooling. It should be possible to reach 
BEC of the dimers and cool the Bose-condensed gas to 
temperatures of the order of its chemical potential. 

Then, converting the molecular BEC into fermionic 
atoms by adiabatically sweeping across the Feshbach res- 
onance to the negative side [29], one obtains the atomic 
Fermi gas at extremely low temperatures T ~ 10~ 2 Tp 
which can be below the BCS transition temperature T c 
[25]. At these temperatures one has a very strong Pauli 
blocking of clastic collisions. The collisional rate is sup- 
pressed as 15(T/Tp) 2 [5], i.e. by a factor of 10 3 , and one 
expects the collisionless regime for the Fermi gas. This 
circumvents the main difficulty for identifying the BCS- 
paired state via a free expansion of the cloud released 
from the trap. This difficulty was met in experiments 
with hydrodynamic Fermi clouds [5] which in cylindrical 
traps are characterized by the same asymmetry of ex- 
pansion below and above T c . Omitting small mean- field 
effects, the collisionless cloud above T c expands symmet- 
rically, whereas in the superfluid state the expansion be- 
comes significantly asymmetric [26]. 
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